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Abstract — This paper deals with the problem of stabilizing 
an unstable system through networked control over the wireless 
medium. In such a situation a remote sensor communicates the 
measurements to the system controller through a noisy channel. 
In particular, in the AWGN scenario, we show that protograph- 
based LDPC convolutional codes achieve anytime reliability and 
we also derive a lower bound to the signal-to-noise ratio required 
to stabilize the system. Moreover, on the Rayleigh-fading channel, 
we show by simulations that resorting to multiple sensors allows 
to achieve a diversity gain. 

I. Introduction 

In the field of control theory, it is of growing interest 
the study of networked control systems, where the mea- 
surement sensor and the controller are not physically co- 
located. Such a case is suitably modelled by supposing 
that the remote sensor transmits its measurements to the 
controller through a noisy communication channel. From 
the information-theoretical point-of-view, such communication 
problem has many differences from the ordinary reliability 
problem on a point-to-point link. Such differences arise es- 
sentially from the fact that systems must be controlled in 
real-time, while the usual approach does not consider delay 
as a primary parameter. Moreover, past decoding errors at 
the receiver may have a catastrophic effect if they are not 
eventually corrected as time proceeds. 

Based on the previous considerations, Sahai and Mitter [ 1 ] 
introduced the new concepts of anytime reliability and of 
anytime capacity. Loosely speaking, an encoding-decoding 
scheme is said to be anytime reliable if its bit error probability 
decreases exponentially with delay d, i.e., is proportional tcfM 
e~P d , where /3 is the anytime exponent of the scheme. Then, 
the anytime capacity C((3) is the supremum of achievable rates 
for schemes with anytime exponent j3. For further information, 
see |2| and references therein. 

Anytime-reliable nonlinear tree codes were first proven 
to exist in |5J and then further developed in [4]. Random 
linear codes were first introduced in [5]. Later, Sukhavasi 
and Hassibi [6| showed that causal random linear codes with 
maximum-likelihood (ML) decoding are anytime reliable with 
high probability. Such schemes are characterized by a high 
decoder complexity, although in J6) a decoder with reasonable 
complexity is proposed for the erasure channel. In [7] Dossel et 
al. proposed a low-density parity-check (LDPC) convolutional 

'in this paper, we use the Euler number c instead of 2 as the base of the 
exponential. 



encoding scheme on the erasure channel which is shown to 
be anytime reliable. In this case, a belief-propagation decoder 
allows for achieving anytime reliability at an affordable com- 
plexity. 

Since the erasure channel is more suitable for modeling 
the behavior of upper layer communications, in this work, 
we study anytime reliable LDPC convolutional codes over the 
wireless channel. The rationale behind this choice is that in this 
case LDPC decoding algorithms can make use of soft infor- 
mation and thus are expected to be more efficient. Moreover, 
a physical-layer encoding-decoding scheme can better exploit 
the potentialities of a wireless multi-node control network 
where a network of remote sensors transmit simultaneously 
their measurements to the system controller. 

The contributions of the paper are as follows: 

• We generalize the design of the LDPC convolutional 
codes proposed by |7j to a broader class. 

• We prove that such class of codes are anytime reliable 
on the AWGN channel, and we give a lower bound on 
their anytime exponent. 

• We show by simulations that, in the presence of fading 
channel between the sensor and the controller, better 
stability margins can be achieved by using multiple 
sensors. 

The paper is organized as follows. In Section[ll] we describe 
the (single-sensor) model of the considered system. In Section 
[III] we generalize the LDPC convolutional codes introduced 
in J7). In Section IV we derive a lower bound on the anytime 



exponent for the LDPC convolutional encoding scheme on the 
AWGN channel. Eventually in Section [VJ through numerical 
analysis, we validate the theoretical results obtained in the 
previous sections and we show simulation results for the sensor 
network scenario. 

II. System model 
We consider the discrete-time dynamic time invariant sys- 
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where x t £ R" 1 is the state of the system at time step t, A 
and B are n x x n x and n x x n u real matrices, respectively, 
u t is the control input, and w t is a zero-mean bounded noise 
process. The system in ([T} is supposed to be unstable, i.e., it is 

2 Column vectors and matrices are denoted by lowercase and uppercase 
bold letters, respectively. 



Dynamic system 



x t+ i = Axt + But + w t 









Controller 




RX 



i r 



Wireless V 1 






Communication L_ 




TX 


Channel I 
I 




V. 





Remote sensor 



yt = Cx f + v t 



Fig. 1. Model of the discrete-time linear dynamic system: remote sensor, wireless communication channel and system controller. 



characterized by p(A) > 1, where p(A) is the spectral radius 
of the matrix A, that is the largest eigenvalue modulus of A. 
The state x t of the linear system in ([T]) is measured by a remote 
sensor (see Figure JTJ providing the size- 

yt = Cx ( 
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where v t G W ly is a zero-mean bounded noise process 
independent of w t . 

The remote sensor is equipped with one antenna and sends 
its measurement to a controller through a noisy wireless 
communication channel. Specifically, in the considered setting 
the sensor at each time step t first quantizes the measurement 
y t into a A: -bit vector q t . In this paper, as in (6), we consider 
a uniform lattice quantizer. The quantized measurements are 
the input of a channel encoder £ t , whose output e t is a binary 
vector of length n. The code rate is thus R c — k/n. The 
encoded bits are modulated to a vector of m symbols, s t , 
belonging to a given constellation and then trasmitted over 
the wireless channel. 

The received signal is given by 



H t s t + z ( 



(3) 



where z f is a size-m vector representing additive noise with 
i.i.d. complex circularly symmetric Gaussian random entries 
with zero mean and unitary variance. The m x m diago- 
nal channel matrix H f = diag(/iti, . . . , ht m ) is such that 
its j-th diagonal element represents the channel coefficient 
experienced at time t by the j-th transmitted symbol, for 
j = l,...,m. 

At time step t the receiver processes the received signal 
r t and obtains the soft estimates of the quantized measure- 
ments. Specifically, r t is first processed by a ML demodulator 
which outputs soft estimates e t of the coded bits e t . The 
soft estimates e t and the estimates computed at time steps 
t = 1, . . . ,t — 1, i.e., ei, . . . , e t _i are then sent to a decoder 
T> t which outputs estimates qolt, • • • , qt|< of the quantized 
measurements. The notation q T | t represents the estimate of 
q r obtained at the receiver at time step t. 

Finally a digital-to-analog converter provides estimates 
yilt) • • • > Ytlt °f me observations. Again the notation y T i t 
represents the estimate of y T obtained at time step t. These 
estimates are sent to the system controller. 



The system controller is in charge of generating suitable 
commands u t in order to keep the system stable. In par- 
ticular at time t the controller takes as input the estimates 
Yilt) • • • > Ytlti an d outputs the control u t . The controller is 
a chain of t filters. The r-th filter produces the output x T | t , 
which is an estimate of the state x T at time step t, and has two 
inputs: the estimate x T _iu and the vector y T \t- The output of 
the t-th filter, x t | t is then providing the estimate of the current 
state. The command u f is finally obtained as a linear feedback 
of the state estimate x t | t 

= Kx t | t 

where the n x x n u matrix K is chosen to stabilize the system, 
i.e. so that p(A + BK) < 1. In this paper, hypercuboidal 
filters [6J have been employed in the system controller. 



III. LDPC CONVOLUTIONAL ENCODING/DECODING 
SCHEMES 

Following [7], we consider a channel coding scheme based 
on systematic LDPC convolutional codes. Precisely, the en- 
coder at time t, £ t , is the systematic encoder corresponding 
to the parity-check matrix 
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where all matrices Zi, i — 0, ... ,i — 1, are (n — k) x n 
sparse binary matrices, and Zo is full-rank (over GF(2)) 
in order for Z[ 1:t j to have full row rank. (Notice that we 
have explicitly restricted our focus to a Toeplitz parity-check 
matrix, for simplicity.) Such coding scheme is causal thanks 
to the lower-triangular structure of the parity-check matrix. 
Moreover, it can be considered as a convolutional code with 
infinite memory, whose number of states is equal to 2(*~ 1 ' fe 
at time t. 

Structurally, as in J?), we have built the parity-check matrix 



Zj 1:t ] starting from a protograph [8] matrix P[i : t], given by error probability of the LDPC encoding/decoding scheme for 
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where all matrices P^, i = 0, . . . , t — 1, are (no — fco) X no 
matrices with nonnegative integer entries and Po is full-rank. 
Zr 1:t ] is obtained by lifting P[i : t] to order r, namely: 

• each zero of Pn. t i is lifted to a r x r all-zero matrix, and 

• each nonzero entry of Pmt] equal to 6 is lifted to the 
modulo-2 sum of 6 permutation matrices of size r x r, 
chosen at random between the r! possible permutation 
matrices of that size. 

In this way, we obtain a (n — k)t x nt sparse binary matrix 
Z[ 1:t ], with (n — fc) = r(no — ko) and n = rno. Notice that 
R c = k/n = fco /no. Since permutation matrices are chosen 
at random, we actually obtain an ensemble of codes, each one 
corresponding to a given choice of the permutations. The fact 
that Z[ 1:4 ] is sparse even if P[i : t] is not, allows us to choose 
the latter with a certain degree of freedom. In particular in our 
work we choose 
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and, for r = 1,2,. 
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where fco = no — fco- Note that the choice of the protograph 
matrices in |6]) and |7]) is more general with respect to the 
choice made in (7J where fco = 1 and no = 2. 

In our system the decoder T> t implements belief propagation 
(BP) [9 1 on the bipartite graph defined by Z[ 1:t ] . In this case the 
advantage of deriving the code by lifting a protograph relies 
on the fact that the local structure of the code graph always 
looks like the protograph one ( JSJ, flO)), while the probability 
of short cycles (detrimental for BP) is reduced by increasing 
r. Thus, convergence properties of BP decoding algorithm can 
be studied directly on Prut], while neglecting the effects of 
cycles. 

As for ordinary LDPC codes |9j, P[i : t] is interpreted as the 
adjacency matrix of the protograph at time t, where columns 
represents variable nodes (VNs) and rows represent check 
nodes (CNs). If a given element of Prut] i s equal to b, there are 
b edges connecting the corresponding CN and VN. Moreover, 
if b > 0, the two nodes are neighbors. Let Af c (i) be the set of 
VNs that are neighbours of z-th CN. Analogously, let Af v (j) 
be the set of CNs that are neighbours of VN j. 

IV. Anytime reliability of LDPC convolutional 

CODES ON THE AWGN WIRELESS CHANNEL 

In this section, we study the anytime reliability of the 
protograph-based LDPC convolutional codes. We first provide 
a theoretical analysis where we derive a bound on the bit 



a generic noisy channel. Then in Section IV-B we specialize 
to the AWGN case. 

A. A lower bound on the anytime exponent of LDPC encod- 
ing/decoding schemes 

Let P e (t,d) be the probability that, at time f, the oldest 
decoding error made by the decoder T> t is d steps back in the 
past: 



P e (t,d) = P{{q t _ d+1 | t + qt-d+i}n 

{4r|t = Qt,t < t-d + 1}}. 



(8) 



We say that the encoding-decoding scheme is anytime reliable 
on a given channel if it satisfies: 



P e (t,d) < Ke~ pd , Vi, d > d 



(9) 



where K, f3 and do are positive constants that depend on the 
coding scheme and on the channel. If the code satisfies (|9j, 
then /3 is called its anytime exponent (on that channel). 

Sukhavasi and Hassibi derive in [6| the conditions under 
which an anytime reliable encoding-decoding scheme can be 
used to stabilize the system of in the mean-square 

sense, so that the expected value of ||x f || 2 is bounded for 
all t. It is shown in [6] that, when using hypercuboidal filters, 
mean-square sense stability is achieved by a code with anytime 
exponent (3 satisfying (3 > 2 log p(A), where A is the n x x n x 
matrix whose elements are the absolute values of the elements 
of A. 

In order to assess the performance of the BP decoder, in 
what follows, we slightly modify the P-EXIT approach of 
jTJ, p0) . In particular, we suppose that the BP messages 
exchanged between VNs and CNs are sent through AWGN 
channels, and we track the evolution of the SNR of such 
channels with the iteration index of the BP algorithm. 

Let us define the following variables: 

• Pch(j)- the physical-channel SNR for VN j. 

• Pq^y t {i, j): the SNR for message travelling from CN i 
to VN j at the l-th iteration of the BP algorithm and at 
time step t (if i G N v (j)). 

• t {i, j)' the SNR for message travelling from VN 
j to CN i at the Z-th iteration of the BP algorithm and at 
time step t (if i G N v (j)). 

Then, the approximate SNR evolution at time t can be de- 
termined through the following set of update equationsp] JlT) : 

• Initialization: For j = 1, . . . , not, i g Af v (j): 
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CN to VN update: For i = 1, . . . , k t, j e Af c (i) 
( \ 
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For simplicity, the update equations are given in the hypothesis that the 
protograph matrix is binary. However, the results hold in general. 



• VN to CN update: For j 
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. Output decision variable SNR: For j = 1, 
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The function M(p) appearing in ( 11 1 is defined as M(p) 
J- 1 (1 - JO)), where 
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gives the mutual information between the input of a binary- 
input AWGN channel with SNR p and the corresponding 
output. Notice that M (p) is a nonnegative, strictly decreasing 
function of p and that M~2_(o) = M(p). 

It can be easily proven Jill that the sequences p^Q^ vt {i,j) 

and Pfj^y t (i, j) are monotonically increasing with iteration 
index I. Moreover, these sequences are bounded, as long as 
Pch(j) i s bounded for all j. Thus, they converge to a limit 
when I goes to infinity. Let us call such limits p^\y t {i-,j) 
and py°lc which are a function of the channel SNR 

values. The output decision variable SNR for VN j, after a 
large number of iterations, is then given by 
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Next, starting from ([8]) we compute an upper bound to 
P e (t, d) after a large number of BP iterations as follows: 
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where (a) follows from ([8j, (b) follows from the fact that 
the code is systematic and the union bound and (c) from the 



Chernoff bound on Q(x). Thus, thanks to (|9j) and (16 1, the 
anytime exponent of the considered coding scheme can be 
lower-bounded by 



P>P= lim 
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Thus, a sufficient condition for the stabilization of system 
(|l|l-(|2ji in the single-node scenario is that (3 > 21ogp(A). The 
value of /? can be obtained numerically, for a given coding 
scheme, thanks to ([T0]i-([T3]l. 



B. The AWGN case 

For the AWGN case, p c h{j) = Pch f° r a H j- Notice that, 
because of (|6|l and (|7j, every systematic VN is connected to a 
CN that is connected to a degree- 1 nonsystematic VN. Since 
the message coming from a degree- 1 VN is set to p c h at 
every iteration, thanks to the fact that M (p) is monotonically 
decreasing with p, we can upper-bound the CN-to-VN message 
exchanged at iteration I as 
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which clearly holds also for I — > oo. By plugging the above 



bound into (15 1, we can write the following upper bound on 



the output SNR for systematic variable j: 

p { r\3)<mu)\ + i)PcH. 



(19) 



It is shown in fTT) that this upper bound is actually reached 
for p c h — > oo if some mild conditions are satisfied. 

Thus, for sufficiently large p c /j, the lower bound on the 
anytime exponent can be approximated by 
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In order to achieve anytime reliability for sufficiently large 
Pch, me VN degrees must then increase linearly with d. If 



\Af v ((t- d)n + =jd + 
P = IPch 



o(d) for all t, then 



(21) 



The above results tells that, in order to stabilize the system in 
([T]i-(|2| over an AWGN channel with the LDPC convolutional 
encoding-decoding scheme and hypercuboidal filters, it is 
sufficient that the channel SNR satisfies p c h > 41ogp(A)/7. 

Notice also that, from the analysis, it seems beneficial to 
increase 7, which corresponds to using a denser protograph: 
however, while this is true for a lifting order r going to infinity, 
increasing 7 may affect negatively the performance for a finite 
r, due to the increased probability of finding short cycles. 



V. Results 

To validate the theoretical results, 
system characterized by 
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with p(A) = 1.997, and with the matrices B and C chosen as 
in (6] Example 1]. Moreover we used the LDPC convolutional 
code of |7j, with ko = 1, no = 2 and P r p = 1 for all r, so 
that 7 = 1. 

If the channel between sensor and controller can be mod- 
elled as AWGN, then, by using (21 1, we conclude that the 



system can be stabilized in the mean-square sense as long 
as Pch > 4 log /o(A) = 4.4 dB. Fig. \2\ shows the error 
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Fig. 2. Anytime exponent of the LDPC convolutional codes on the AWGN 
channel: theory versus simulations. 



Fig. 3. Performance of multiple sensors in the fading scenario. 



probability P e (t, d) versus the delay d. The experimental 
curves (solid lines) obtained by simulations, are compared with 
the theoretical ones (dashed), obtained according to ( p~6] >, in the 
hypothesis that the upper bound of (19i is actually achieved. 
As it can be seen, for p c ^ > 4 dB, the slope of the simulated 
P e {t, d), which corresponds to the anytime exponent, is larger 
than the theoretical one, as predicted from the analysis. Monte 
Carlo simulations actually show that the system is controlled 
for p ch = 4.5 dB. 

A. Extension to the fading scenario 

Through simulations, we have investigated an extension of 
the control system of Section [TT] to the fading scenario in the 
case of a network made of multiple sensors. More precisely, 
we consider the case where there are N identical remote sen- 
sors, whose measurements are subject to independent bounded 
noise. At time t, the i-th sensor, i = 1, . . . , N, obtains the 
size-rij, measurement 
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% = 1,..,,N, are zero-mean bounded noise 
processes independent of each other and of w t . The i-th 
sensor encodes the information as described in Section Hill and 

(i) 

transmits the symbol vector sj to the common receiver, which 
is equipped with N antennas. The signal received at the j-th 
receive antenna, j = 1, . . . , N, will then by given by 
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(24) 



where the diagonal channel matrix contains on its diag- 

onal the channel coefficients from sensor i and receive antenna 
j. We assume Rayleigh fading, so that the instantaneous SNR 
at the receiver is a random variable exponentially distributed 

as 

f P M= (25) 

Pch 



The receiver performs jointly optimal demodulation of the 
N superimposed transmitted signals. After demodulation, N 
decoders work in parallel to decode the information sent by 
the N sensors. The controller computes the feedback signal u t 
by putting together the reconstructed measurements from all 
sensors. In Fig. [3] we show the performance of the simulated 
control network. We have measured the probability, pioa, 
that, after 100 time steps, the Euclidean distance between the 
system state and the estimated state is larger than 10 3 , as a 
function of the average SNR, ~p ch , for N = 1,2, 3. The power 
transmitted by each sensor has been normalized so that the 
total transmitted power is the same for the three cases. We have 
used a lifting order r = 12 in the LDPC code construction. 
As it can be seen, pioo decreases with p ch . Notice also that 
the slope of pioo increases with N. We can conclude that a 
diversity gain can be obtained by using multiple sensors. 
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